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Evolution of a network of vortex loops in Hell due to the fusion and breakdown of vortex loops is 
studied. We perform investigation on the base of the "rate equation" for the distribution function 
n(l) of number of loops of length I proposed by Copeland with coauthors. By using the special 
ansatz in the "collision" integral we have found the exact power-like solution of "kinetic equation" 
in stationary case. That solution is the famous equilibrium distribution n(l) oc Z~ 5//2 obtained 
earlier in numerical calculations. Our result, however, is not equilibrium, but on the contrary, it 
describes the state with two mutual fluxes of the length (or energy) in space of the vortex loop sizes. 
Analyzing this solution we drew several results on the structure and dynamics of the vortex tangle 
in the superfluid turbulent helium. In particular, we obtained that the mean radius of the curvature 
is of order of interline space. We also obtain that the decay of the vortex tangle obeys the Vinen 
equation, obtained earlier phenomenologically. We evaluate also the full rate of reconnection events. 

o ' 1 

A network of one-dimensional singularities appears in various physical systems affecting important properties of 
them. As examples we would point out quantized vortices in quantum fluids both in turbulent regime (see e.g. 
book [[J and papers Q, 0L 0) an d thermodynamically equilbrium state (J3, 01). Other examples are the flux 
tubes in superconductors [g, dislocations in solids 8j, global cosmic strings polymer chains |9g . To evaluate 

contribution of mentioned structure in various effects, such as thermodynamic and kinetic properties, phase transition 
etc. we have to know structure and dynamics of network of singularities. 

The whole evolution of network of chaotic set of line consists of two main ingredients. The first is the motion of the 
individual elements of each of lines, due to specific equations of motion, different for each of cases listed above. For 
instance elements of vortex filaments in turbulent superfluid 4 He move obeying the Biot-Savart law supplemented 
by the friction force and the external fiow/counterflow if any. Besides of motion of elements there is another very 
1 ^ 1 ■ important element of the whole dynamics, common for all systems - fusion and breakdown of the loops or recombination 
of open strings due to reconnection processes. Further for dcfinitcncss we will talk about vortex loops in superfluid 
■ helium. 

It is widely appreciated that the reconnection processes influence both the structure and evolution of vortex tangle. 
For instance, Feynman in his pioneering paper devoted to superfluid turbulence proposed scenario how the vortex 
(~~-» tangle decays in absence of the external counterflow. According to this scenario a fusion of small vortex rings into 
larger ones as well as a breakdown into smaller ones is possible at the moments of the reconnection events. In 
assumption that on the average the last property dominates, the cascade like process of formation of smaller and 
smaller loops forms. When the scale of the small rings becomes of the order of the interatomic distances, which is 
the final stage of the cascade, the vortex motion is degenerated into thermal excitations. The idea that degeneration 
of the vortex tangle occurs due to cascade-like transferring of the length in space of scale of sizes of vortex loops was 
indirectly confirmed only in numerical calculations, where the procedure of artificial elimination of small loops had 
1 ■ been used 

In spite of the recognized importance of the reconnecting loops kinetics, the numerical results remain main source 
of information about this process. The obvious lack of theoretical investigations interferes with the deep insight 
in the nature of this phenomena (this situation had been recently discussed in |l7|). For instance it is not clear 
how the cascade of length in space of vortex loops sizes is formed, what mechanisms are responsible for this, what 
. £^ ] quantities determine intensity of cascade, and why at all the breakdown of the loop prevail. The scarcity of analytic 
investigations related to incredible complexity of the problem. Indeed we have to deal with a set of objects with 
not fixed number of elements, they can born and die. Thus, some analog of the secondary quantization method is 
required with the difference that objects (vortex loops) themselves possess an infinite number of degree of freedom 
with very involved dynamics. Clearly this problem can hardly be resolved in nearest future. Recently in 10] much 
more modest approach, based on the "rate equation" for distribution function n(l) was elaborated in context of cosmic 
strings. Following this work we introduce distribution function n(/,t) of density of loop in "space" of their lengths. 
It is defined as the number of loops (per unit volume) with lengths lying between I and I + dl. Due to reconnection 
processes n(l,t) can vary. We discriminate two types of processes, namely the fusion of two loops into the larger 
single loop and the breakdown of single loop into two daughter loops. The kinetic of vortex tangle is affected by 
the intensity of the introduced processes. The intensity of the first process is characterized by the rate of collision 
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A(li, l 2 , 1) of two loops with lengths l\ and l 2 and forming the loop of length I = lx + l 2 . The intensity of the second 
process is characterized by the rate of self-intersection B(l,li,l 2 ) of the loop with length I into two daughter loops 
with lengths l x and l 2 . In view of exposed above we can directly write out the master "kinetic" equation for rate of 
change the distribution function n(l,t). 



dn(l,t) 



A{l u l 2 ,l)n{h)n{l 2 )8{l-h-l 2 )dhdl 2 h + h^l (1) 



at 

' A(h,l,l 2 ,)S(l 2 -h-l)n(l)n(h)dhdl 2 k + l^l 2 (2) (I) 

A{l 2 ,l,l u )Kh-h-l)n{l)n{h)dhdh h + l^h (3) 

B{l x ,l 2 ,l)n{l)5{l-h-l2)dhdl 2 l^k + l 2 (4) 



BQj^hWh-l-hMl^dhdh h^l + l 2 (5) 

BUJxJ^Sih-l-hM^dhdh . l 2 ^l + h (6) 

All processes are depicted at the left of each line. Coefficient A and B were evaluated on the base of qualitative 
considerations in |10| . A bit more rigorous way, having universal character was elaborated in work |18| . Both models 
give similar expressions of type 

A(h > h,l)=b m V l l 1 l 2 , B(h,l-h,l) = b s ™ (2) 

Here Vi is some characteristic velocity ov elements of line, b m and b s are some constants depending on the model. 
For instance in work on evolution of network of cosmic strings autors offered values about 0.1 -j- 0.3 for both 
coefficients (not necessary equal) . In paper ^ij on vortex loops in superfluid turbulent Hell there was offered b m ^ 1/3, 
b s ~ 0.0164772. The quantity £0 appears in as the persistency length of theory of randomly walking chains modelling 
chaotic cosmic st ring s. For the vortex tangle in superfluid helium this quantity is associated with the mean radius 
of curvature (see |l9(). Both these approaches are failed for scales near £ j therefore usually this value appears as a 
low cut-off. Equation Q had been studied (mainly numerically) in papers (lo| , [llT | , where some conclusions about 
evolution of network of cosmic string were made. 

In the present work we demonstrate that the master "rate equation" Q has exact stationary power-like solution of 
form n(l) = C *l s . We discuss the physical meaning of this solution and apply it to describe some properties of vortex 
tangle in the superfluid turbulent helium. To find power-like solution of form n(l) = C*l s we use the Zakharov ansats, 
which is the special treatment of the " collision" integral in equation (QJ . This trick was elaborated by Zakharov for 
the wave turbulence (see e.g. |2fij|). now we will show how it works in our case. Let us take for instance first and 
second integrals in the "collision term" of (Q. Let us further perform in the second integral the following change of 
variables. 

« = ,,(<), ^r.(i), ,=,(<). (3, 

Under this change of variables various factors in the integrand of the second integral transforms as follows 

r- 1 



6(h-h-l)-> [■=-) 5{l-h-l 2 ). 
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As result the second integral in the "collision" term takes a form (additional term 3 in the power counting appears 
from the Jacobian of transformation) 

, \ 2+2s-l+3 

A{h~l 2 , l)n(h)n(l 2 )5(l - h - h)dhdl 2 . (4) 

It is easy to see that the transformed second term in the "collision" integral in rhs of master kinetic equation Q 

turns into first integral with additional factor I l/l 2 ) in the integrand. Performing the same procedure for all lines 
we conclude that the "collision integral" of the "rate equation" can be written as 

^ \ 4+2s / ^ \ 4+2s^ 



A(l u l 2 ,l)n(h)n(h) yi- y-^J J S(l - h - h)dhdl 2 (5) 

^ v s+3/2 , , x s+3/2 N 



B(l u l 2 ,l)n(l) \1- {j-j j 5(1 - h - h)dhdl 2 

/ , \ 4+2s , ,\4+2s / s s+3/2 , x s+3/2 

For s — —5/2 both expressions 1 — f ^ J — (^-J andl— f ^ ) — (ll) are equal to (^ — ^i —^ 2 )/^- Thus 
the integrands of both integrals in 10 include expressions of type (x) S(x) and these integrals vanish. This implies in 
stationary case the power-like solution n = C * l~ b / 2 for distribution function n(l,t) of density of loop in "space" of 
their lengths takes place. 

Let us discuss the physical meaning of the solution obtained. First of all we stress that it does not relate to detailed 
balance, it rather corresponds to nonequilibrium state. To clarify what is happening we introduce density L(t) (in 
space of sizes I) of full length per unit of volume 

. j , , . total length 

L(l,t) =l*n(l,t) = — t — — -. (6 

unit 01 volume^mterval 01 length 

The total length (per unit volume), or the vortex line density 

C(t) = [ L(l,t)dl= [ l*n(l,t)dl 



is obviously conserved during the reconnections events dC(t)/dt = 0. Conservation of the vortex line density can be 
expressed in form of continuity equation: 

This form of equation states that the rate of change of length is associated with " flux" of length in space of sizes of 
the loops. Term "flux" here means just redistribution of length among the loops due to reconnections. Expression 
for P(l) is obtained by multiplication kinetic equation by I and by rewriting the "collision" term in the shape of 
derivative with respect to I. Result is (substitutions — x and l 2 /l — y had been used below) 

P = (P^A I I \b m V lX yC 2 x s y s f 1 - ( -) - (-) ) 6(1 - x - y)dxdy (8) 



5 + 2s J J J 2 y \ \x_ 

is+5/2 \ r r\ 1 / /1 \ s + 3 /2 / 1 \ s + 3 / 2 \ 

1 " ^sVlj^^C 1- - - - )S(l-x-y)dxdy. 



5 + 5/2/7 7 2 s (fox) 3 / 2 I \xj \y 

Both integrals in relation © coincide with integrals in therefore they vanish for s = —5/2. However they have 
pre-integral factor with the denominators, which also vanish for s = —5/2. Calculating indeterminacy 0/0 we obtain 
the final expressions for "flux" of length in space of sizes of the loops 

D 12.555^ 2 , 5.545 

P = ~^—C 2 b m Vi - ^pCbeVx (9) 
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The positive sign of the first corresponds to flux of length in direction of large scales. This is justified, since the 
fusion processes lead to formation of large loops. The negative sign of the second term corresponds to flux of length 
in direction of small scales. This is justified, since the breaking down processes lead to formation of small loops. 
The approach elaborated above allows to draw several conclusions concerning the structure and dynamics of real 
vortex tangle in turbulent Hell. To do it we have to specify quantity V] , which enters into the rates coefficients of 
the both merging and breaking down processes. Since the averaged radius of curvature is £o, we estimate the velocity 
factor Vi to be of order of «/£o (k is the quantum of circulation). Thus the only parameters of the whole theory are 
the quantum of circulation k and the mean radius of curvature £o ■ 

Vortex Line Density. In steady case the positive flux of length in exactly compensates the negative flux. 
Equating these two terms we have 

5.455 b s 1 1 

c - u^k:^- Cvld ^- (10) 

New numerical parameter Cvld ~ 0.2 I . Thus the power -like solution of the master Kinetic equation n(l) is 

, n Cvld ,-5/2 

n (0 = -37a- i 1 ( n ) 
so 

Accordingly the total length L per unit of volume (we recall that quantity £o serves as the low cut-off ). 



c=l /*»(/.,// = |^p. n->) 



Result <|I2[1 is remarkable. The idea that interline space S = C x ' 2 is of order of mean radius of curvature £o was 
launched by Schwarz Earlier it was confirmed only in numerical simulations and the nature of this phenomenon 
was not clear. We proved that this relation appears due to kinetics of colliding vortex loops. 
Decay of Vortex Tangle (Vinen Equation). Let us suppose that there is some sink for loops of small sizes. 
Then the steady situation is reached when there is some source which generates length. This is the mutual friction 
for the superfluid turbulence. If we switch off the source of length, the total length will decrease. The attenuation 
of vortex line density is related to negative flux P neg (in direction of small scales) of the full flux and can be 
evaluated from the continuity equation for density of length L(l, t) in space of the loop sizes, which we integrate over I 

dL(l,t) dP{l) _ dC(t) 
dt + dl ~ dt 

Evaluating P neg from relation for full flux Jjj), where constant C is used from l|10(l we conclude that 

^ — ~Cvinen K £ 2 ■ (13) 

Here Cvinen is new numerical factor. For parameters, which we adopt here Cvinen ~ 0.05 -r 0.2. Expression l |13| l is 

the famous Vinen equation [3§ obtained him from experiment. It is remarkable that this relation appears due to the 
reconnection processes. The own dynamics of filament specific for various systems is absorbed by the persistency 
length £o , which has dropped out the Vinen equation at all. Thus the (|I 3|l has universal character and can be applied 
for other systems, e.g. for cosmic strings. We also would like to note that our calculations confirmed the brilliant 
Feynman's conjecture on the formation of cascade-like breakdown of vortex loops. 

The Full Rate of Reconnection. The full rate of reconnection N rec can be evaluated directly from master "rate 
equation" Indeed, this equation describes change of nil) due to reconnection events. It takes into account sign 
of events, depending on whether the loop of size I appears or dies in result of reconnection. Therefore, if we take all 
terms in collision integral with the plus sign we obtain the total number of reconnections. The according calculations 
lead to result 



rec 



1 K.{b s Cy LD + b^C'vLD) _ n r5 / 2 

d 4"o 



where C rec one more constant of order 0.1 — 0.5. That results agrees with the recent numerical investigation by 
Barenghi and Samuels 
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The problem of nonequilibrium dynamics of vortex loops, which merge and break down due to reconnections has been 
considered. The description was performed on the base of the "rate equation" for distribution of number of loops in 
space of their lengths. By use of special substitution of variables in the collision integral (Zakharov ansatz) we found 
the power-like solution of the kinetic equation. That is non-equilibrium solution characterized by two mutual fluxes 
of length (energy) in space of loop sizes. 

The result obtained were used to draw some conclusion about the structure of the vortex tangle appeared in the 
superfluid turbulent HcII.(i) In particular we obtained that vortex line density is of order of inverse squared mean 
curvature of lines, (ii) We also found that the rate of decay of the vortex line density is proportional to squared 
vortex line density itself multiplied by quantum of circulation and numerical factor of order of unity, (iii) we also 
estimate the total number of reconnections, which turned out to be of order of quantum of circulation multiplied by 
the vortex line density in power 5/2. 

I am grateful to participants of the workshop "Superfluidity under Rotation" (Manchester, 2005) for useful discussion 
of the results exposed above. This work was partially supported by grant N 03-02-16179 of the Russian Foundation 
of Fundamental Research. 
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